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XMM-Newton observation of the eclipsing binary Algol 
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ABSTRACT 

We present an XMM-Newton observation of the eclipsing binary Algol which 
contains an X-ray dark B8V primary and an X-ray bright K2IV secondary. The 
observation covered the optical secondary eclipse and captured an X-ray flare that 
was eclipsed by the B star. The EPIC and RGS spectra of Algol in its quiescent 
state are described by a two-temperature plasma model. The cool component has 
a temperature around 6.4x10^ K while that of the hot component ranges from 2 
to 4.0x10'^ K. Coronal abundances of C, N, O, Ne, Mg, Si and Fe were obtained 
for each component for both the quiescent and the flare phases, with generally 
upper limits for S and Ar, and C, N, and O for the hot component. F-tests 
show that the abundances need not to be different between the cool and the hot 
component and between the quiescent and the flare phase with the exception of 
Fe. Whereas the Fe abundance of the cool component remains constant at ~0.14, 
the hot component shows an Fe abundance of ~0.28, which increases to ~0.44 
during the flare. This increase is expected from the chromospheric evaporation 
model. The absorbing column density Nh of the quiescent emission is 2.5x10^° 
cm~^, while that of the flare-only emission is lower than 1.3x10^° cm~^. Together 
with previous observations, this Nh difference suggests that there exists a belt 
of cool gas above the corona and the belt runs parallel to the equator of the K 
star. 

Subject headings: stars: abundance - star: binaries: eclipsing - stars: flare - 
stars: individual: Algol 
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1. Introduction 

Algol {(3 Per) is a nearby eclipsing binary with an early main- sequence (B8V) primary 
and a cool subgiant (K2IV) secondary. This system is only 28.46 pc away (Hipparcos parallax 
measurements, ESA 1997) and has a period of about 2.87 d and an orbit inclination of 81°. 
The radii of the two companion stars {Rb and Rk) are 3.0 and 3.4 Rq, respectively, and their 
separation is 14.14 Rq (Hill et al. 1971; Richards 1993). Algol was first detected in X-rays 
by SAS 3 (Schnopper et al. 1976) and confirmed as a strong X-ray emitter by Harnden et 
al. (1977). It is generally accepted that the K star has an active corona and accounts for 
most of the X-ray flux of the system, whereas the B star is X-ray dark (White et al. 1980). 
This was proved for the first time by Chung et al. (2004) by the detection of Doppler shifts 
of the spectra caused by orbital motion of the K star. Given an active corona, flares appear 
typically every day with a duration of several to many hours. 

Algol shows elemental abundances which differ from those of the sun (Antunes et al. 
1994; Schmitt & Ness 2004). Antunes et al. (1994) quote abundances of Fe, O, Mg, Si, 
S, Ar and Ca which are lower than the solar photospheric values by a factor of 2-3, and N 
to be less than 0.1. With the exception of N these results have been generally conflrmed 
by Chandra observations (Schmitt & Ness 2004). Schmitt & Ness (2002) studied a sample 
of late type stars including Algol and they found an enhancement of N. Drake (2003) also 
suggested that the N abundance is enhanced by a factor of 3, and C is depleted by a factor 
of 10 (both relative to HR1009, whose C/N abundance is consistent with that of the solar 
photosphere). These C and N abundances would lead to the conclusion that the K star has 
lost at least half of its initial mass. 

That flares have different elemental abundances has been concluded from previous mis- 
sions, such as GINGA (Stern et al. 1992), ROSAT (Ottmann & Schmitt 1996), and Bep- 
poSAX (Favata & Schmitt 1999). They suggest that the elemental abundances of Algol's 
flaring region resemble more the photospheric than the quiescent coronal abundances, which 
might be a consequence of chromospheric evaporation (Giidcl et al. 1999; Favata & Micela 
2003 and references therein). At the beginning of the event, fresh chromospheric material is 
evaporated in the flaring loop(s), and enhances the corona elemental abundances; once the 
material has been transported to the flaring corona structure, the fractionation mechanism, 
which is responsible for the lower abundance, would start operating, bringing the corona 
abundance back to its quiescent value. 

The absorbing column density {Nh) inferred from X-ray observations is usually higher 
than the column density of the interstellar medium (ISM) between Algol and the observer. 
Welsh et al. (1990) gave an upper hmit of the ISM column density of 2.5x10^^ cm~^ using 
the ISM Nal D fine absorption of the B8 primary. Stern et al. (1995) got a similar result 
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from the He/H ratio using EUVE observations. The Nh derived from the ROSAT PSPC 

observation of Algol is about 5x10^^ cm~^ (Ottmann 1994). Favata & Schmitt (1999), using 
BeppoSAX, reported that Nfj increased up to > 10^^ cm~^ during the flare rise and decreased 
during the decay, which may possibly be associated with moving, cool material in the line 
of sight, i.e. a major coronal mass ejection. 

In this paper we present a time resolved spectroscopy of Algol, using XMA'I- Newton 
European Photon Imaging Camera (EPIC) and Reflection Grating Spectrometer (RGS) 
observations of Algol. Compared with the instruments used in earlier studies the EPIC 
cameras are the flrst ones to have a combination of moderate spectral resolution, wide energy 
coverage, and high sensitivity, while RGS provides us with unparalleled sensitivity for high 
resolution spectroscopy in the soft X-ray band. These advantages permit us to perform a 
more detailed diagnosis of the plasma properties and their evolution as far as the corona of 
the K type secondary is concerned. The XMM-Newton observation covered the secondary 
eclipse and an X-ray flare was detected. Schmitt et al. (2003) presented a detailed study of 
the geometry of the flare, and we analyse the X-ray spectra of Algol in both the quiescent 
and flaring states. We can therefore study the mechanism of the X-ray emission of Algol 
in both the quiescent and the flaring states. In § 2, we describe the observation and data 
reduction. The spectral properties of the overall and flare-only emissions are given in § 3. 
The proposed local cool absorbing gas structures are discussed in § 4 and a summary is 
provided in § 5. 



2. Observation and Data Reduction 

Algol was observed with XMM-Newton on 2002-02-12 from 04:02:39 to 18:45:09 UTC. 
We use data from the EPIC-PN (Striider et al. 2001), the (two) EPIC-MOS (Turner et al. 
2001) detectors and the (two) RGS spectrometers (den Herder et al. 2001). The EPIC-PN 
was turned on at 04:42:18 and turned off at 18:34:45 UTC, with a net exposure time of about 
50 ks, while the two EPIC-MOS and the two RGS detectors had a bit longer exposure, from 
04:09:01 to 18:39:28 UTC and from 04:02:39 to 18:42:59, respectively. The EPIC-PN has an 
energy coverage of 0.15-15 keV, energy resolution of 80 eV at 1 keV, and a time resolution of 
73 ms in the full-frame imaging mode in which the detector was operated for this observation. 
The EPIC-MOS has a narrower energy coverage of 0.15-12 keV, a better energy resolution 
of 70 eV at 1 keV, and a time resolution of 2.6 s in the standard full-frame mode that was 
used in this observation. The RGS allows much better energy resolution measurements (i.e. 
2.9 eV at 1 keV) in the soft X-ray range (0.3 to 2.1 keV). Algol was positioned at the aim 
point, and thick filters were used in order to reduce the effect of strong optical loading of 
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the CCD chips. 

The X-ray data were analyzed using the XMMSAS software package. We found that 
there is no significant background flare and thus no time interval except the observation 
gaps had to be ehminated from the dataset. As Algol is very bright, the center part of the 
source has been affected seriously by pile-up for the EPIC data, but not for the RGS data. 
Therefore, in the EPIC data analysis, we used the PN counts extracted from an anmilus of 
21.75 arcsec inner radius and 125 arcsec outer radius. For the MOS inner and outer radii 
of 21.25 and 125 arcsec were taken for the selection annulus, respectively. According to the 
integrated PSFs of EPIC-PN and EPIC-MOS at 1.5 kcV (which is near the peak of Algol's 
emission), photons in the annuli chosen make up about 10% of the total source photons for 
the EPIC-PN and EPIC-MOS detectors (Aschenbach et al. 2000). This factor was taken 
into account when we calculated the fluxes and emission measures. 

The corresponding EPIC PN and MOSl 0.3 to 10 keV and the RGSl 0.3 to 2.1 keV 
background corrected light curves for the selected annuli and the entire observation are 
shown in Fig. 1 with time bins of 100 s each. Using the ephemeris Primary minimum — 
HJD2, 445, 739.0030+2.'^8673285£; (Al-Naimiy et al. 1985; Kim et al. 1989), the observation 
falls in the phase range 0.36-0.56 which includes the optical secondary eclipse. 

In order to study the temporal evolution of the spectrum, we divided the observation into 
6 segments and created a spectrum for each. The integration time is 5000 to 10000 s for these 
spectra. For the EPIC-PN, we only used the single events (PATTERN=0) for the spectral 
analysis, and the background was extracted from an annulus of 490 and 635 arcsec (the area 
contaminated by out-of-time events was excluded). For the two EPIC-MOS, we used both 
single and double events, and the background spectrum was extracted from an annulus of 
490 and 660 arcsec. For the two RGS only the first order data were used. Each spectrum 
was fitted with a two-temperature thermal plasma model (VMEKAL model; Mewe et al. 
1995) in XSPEC. The free parameters are temperatures (Tl, T2), emission measures (EMI, 
EM2) and abundances of C, N, O, Ne, Mg, Si, S, Ar and Fe for each of the two components. 
The abundances are in units of solar abundances given by Anders & Grevesse (1989). The 
WABS model (Morrison & McCammon 1983) was used to take care of interstellar photo- 
electric absorption. In the joint fit of the EPIC and RGS spectra, the RGSl, RGS2 and the 
EPIC emission measures were allowed to be different, because the EPIC counts were selected 
from a restricted annulus as described above, and the corresponding emission measures were 
of course much lower than those of the two RGS. 

To study the properties of the flaring plasma, the flare period was cut in two sections, 
i.e. sections 7^ 4 & 5 (c.f. Fig. 1). We chose the quiescent spectrum of the flrst 8000 s of 
the observation (c.f. Fig. 1) as background spectrum, which is justified by the fact that the 



- 5 - 



X-ray eclipse often starts around phase 0.40 (e.g., Schmitt & Favata 1999), corresponding 
to about 10000 s after exposure start of the present observation. Given that the quiescent 
emission dominated in the RGS energy band, use of the RGS data could not be made because 
of the fairly low signal to noise ratio. Taking into account that the quiescent emission was 
also eclipsed during the flare, we scaled the background spectrum by a factor of 0.9, which is 
the eclipsed fraction of the quiescent emission (Schmitt & Favata 1999). It should be noted 
that we assume that the quiescent emission arises from an isotropic corona of the K star, 
and thus no rotation modulation was taken into account when we dealt with the background 
emission of the flare-only spectrum. The spectral analysis was performed with XSPEC. A 
one-temperature MEKAL model in combination with the WABS model was used. 



3. Spectra 

3.1. The spectra of the overall X-ray emission 

The best-fit model parameters of the 6 spectra as well as their errors are listed in Table 
1 & 2 and plotted in Fig. 2, 3 & 4. The errors are calculated in XSPEC taking into account 
the full error projection, but they represent just the statistical errors at 90% confidence. As 
an example. Fig. 5 shows the model fit to the spectrum #4. 

From Tables 1 & 2 and Fig. 2, 3 & 4, we derive the following results: 

(I) The best-fit Nh values are about ~ 2.5 x 10^° cm~^, and remain nearly constant 
during the entire observation. This Nh is a factor of 100 higher than the upper limit of 
the ISM HI + H2 column density toward Algol derived by Welsh et al. (1990). Fitting the 
spectra with A'^^ fixed to 2.5x10^^ cm~^ ends up with higher values. With these two 
sets of Nh, ^^'i degrees of freedom, we can carry out the F-tcst using the ftest-routine 
in XSPEC. The results show that most of the spectra have of a probability of <0.01% with 
Nh that small. An Nh value increase during the fiare rise, which has been detected with 
BeppoSAX (Favata & Schmitt 1999) has not been detected in our observation, i.e., there 
is no evidence for major mass ejection during the flare to the extent that extra absorption 
occurs. We also tried to flt each of the two spectral components with different Nh values 
but the flt did not improve signiflcantly. i.e., the Nh values of the two components show no 
signiflcant difference. 

(II) Both the temperature and the flux of the hot component show signiflcant temporal 
variations similar to the fight curve. This conflrms the ROSAT PSPC results that the increase 
of the count rate during the flare is predominantly the result of an emission measure increase 
of the hot component (Ottmann 1994; Ottmann & Schmitt 1996). 
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(III) Table 3 demonstrates that the abundance pattern is generally compatible with 
that obtained from ASCA (Antunes et al. 1994) and Chandra observations (Schmitt & Ness 
2004). An exception is observed for oxygen abundance of the hot component. But oxygen, 
as well as carbon and nitrogen should be fully ionized at these high temperatures so that 
emission lines cannot reliably be used to constrain abundances, which means that in our 
case the abundances are poorly determined. The relatively high N abundance and N/C 
ratio (for the cool component), which has been detected earher by Schmitt & Ness (2002), 
Drake (2003) and Schmitt & Ness (2004), are clearly present in our observation as well. They 
therefore support the suggestion made in the papers quoted above that the N enhancement 
is related to the CNO-cycle operating in the core of the star. Furthermore, the active corona 
of Algol seems to be Fe-poor, which has been suggested before by Antunes et al. (1994) and 
Stern et al. (1995). Iron is likely to be depleted in the corona (Drake 2003). 

Comparing the abundances of the cool and the hot component (c.f. Table 2, Fig. 3 
& Fig. 4) we note that apart from perhaps the light elements C, N, O there appears to 
be no significant differences for Ne, Mg and Si. Because of the large uncertainties also the 
abundances of S and Ar are not really accessible, but upper limits can be provided. A clear 
difference, however, exists for Fe (c.f. Fig. 6). In fact, the F-test shows that we cannot 
rule out that the abundances of the cool and the hot component are identical but with the 
exception of Fe. Fig. 6 shows the best-fit and the corresponding error contours for the Fe 
abundances of the cool and hot component of spectrum #1, i.e. during the quiescent phase. 
Up to the 4(7 level the abundances differ from each other. For the best fit the abundances 
are Fe(cool)~ 0.14 and Fe(hot)~ 0.28. This result does not only hold for spectrum #1 
but for the entire quiescent phase including the post-flare phase covered by spectrum ^6. 
From the ROSAT PSPC observation (Ottmann 1994) it was concluded that the low and 
high temperature components are rather co-spatial with respect to the stellar surface. We 
therefore speculate that the cool and hot components might come from different layers or 
that the hot regions are much smaller in size than the radius of the K star. We admit that the 
two components are not likely to represent physically distinct temperatures. A temperature 
gradient in the corona might exist, or the two temperatures are just a short-cut expression 
of a temperature dependent emission measure. 

(IV) The abundance of Fe and possibly those of Ne and Si of the hot component appear 
to change during the flare (c.f. Table 2). The Fe abundance increased from 0.25 times 
the solar value in the quiescent state to about 0.5 during the flare. To put this in a more 
quantitative context we carried out an F-test and assess whether the elemental abundances 
are identical for the quiescent and the flare-only emission. This hypothesis can indeed be 
accepted but only if the abundances of Fe between the hot components differ. The cool 
component is consistent with one set of Fe abundance for the quiescent and the flare-only 



- 7- 



phase. Fig. 7 shows the best fit and the error contours for the Fe abundance of the hot 
versus the cool component for a flare spectrum. The Fc-abundances are inconsistent with 
each other well beyond the 5cr level. Whereas the Fe abundance of the cool component stays 
at ~0.14 the best fit for the Fe abundance of the hot component is now at ~0.44, which 
means an increase by a factor of ~1.6 between the quiescent and the flare phase. This Fe 
abundance increase during a flare may be attributed to chromospheric evaporation (Favata 
& Micela 2003). The abundances of Ne and Si for the hot component may show a similar 
behaviour although with much larger uncertainty. This could further support the idea of 
chromospheric evaporation. The elemental abundances of the cool component remain nearly 
unchanged during the entire observation, which is additional support for the hot component 
dominating the flare. 



3.2. The flare-only spectra 

The spectra of the flare-only emission can be well fltted with a one-temperature MEKAL 
model. Table 4 gives the fltted parameters and Fig. 8 shows the model flt to the flare- 
only spectrum #2 (corresponding to the overall spectrum #5 ). The temperature of the 
flare emission is around 6.2x10^ K and it shows no detectable variations during the flare. 
Obviously the flare temperature is much higher than the temperature of the hot component of 
the quiescent state. The elemental abundances are comparable to those of the hot component 
and they show no variation either (see Table 4). The best-flt Nh of the flare emission is 
consistent with no absorption at all for each of the two spectra. Taking into account the 
uncertainties it is still signiflcantly lower than the Nh values associated with the quiescent 
emission, which is consistent with the results of the ROSAT observation (Ottmann 1994; 
Ottmann & Schmitt 1996). We then used the F-test to examine whether the Nh difference 
between the flare emission and the quiescent emission is real or just a statistical fluctuation. 
This test is based on the comparison of the values and the degrees of freedom of the model 
fits to the flare-only spectra with Nh as a free parameter against Nh being fixed to the mean 
Nh of the quiescent emission (2.5 xlO^'^ cm~^). As shown by the values listed in Table 4, 
the Nh of the flare emission is lower than that of the quiescent emission at a signiflcance 
level greater than 99.5%. A similar behavior is suggested by the ROSAT observation with 
an overall Nh of about 5x10^^ cm~^ (Ottmann 1994), while the flare-only Nh is around 
1x10^^ cm-2 (Ottmann & Schmitt 1996). 
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4. Cool absorbing gas surrounding the K star 

The mean Nh value of the quiescent emission derived in this paper is about 2.5x10^'^ 
cm~^, which is at least two orders of magnitude higher than the interstellar HI + H2 column 
density reported by Welsh et al. (1990). The low energy calibration of XMM may affect our 
results, but not as much as two orders of magnitude. Stern et al. (1995) estimated the Nh 
to be 2x10^^ cm~^ using EUVE observation, which is consistent with that of Welsh et al. 
(1990). The A^^ of 5x10^^ cm-^ obtained in the ROSAT PSPC observation is lower than 
our result (Ottmann 1994), but still 20 times higher than the ISM column density. However, 
since Welsh et al. (1990) inferred the ISM column density using the Na line absorption by the 
continuum of the B type primary that is X-ray dark, and since the values obtained in X- 
ray observations correspond to the absorbing column densities of the X-ray emitting K type 
secondary, the higher Nh values thus imply that there exists cool gas just surrounding the K 
star. Furthermore, the result of Stern et al. (1995) was based on the He/H ratio and with no 
ionized He. This would be the case for our model of ionized local absorbing material, which 
would make the EUV transparent but not to X-rays. Because no orbital (and so rotational) 
modulation of Nh has been detected by either ROSAT (Ottmann 1994) or XMM-Newton, 
the cool gas should be distributed quite homogeneously relative to the longitude of the K 
star surface. The A'^^ difference between the ROSAT and XMM-Newton observations might 
further suggest that the amount of cool gas surrounding the K star changes over a long term. 

We consider two possible patterns of the distribution of the cool gas around the K 
star. One pattern consists of a layer which is co-spatial with the corona. As shown in Fig. 
9(a), the very low Nh of the flare-only emission means that the bulk of the flare does not 
suffer from absorption by the cool surrounding gas, and therefore the height of the cool gas 
is smaller or at most comparable with the height of the flare, i.e. ~ O.IRk cm (Schmitt 
et al. 2003). Taking into account that the height of the corona is comparable with the 
stellar radius (White et al. 1986; Ottmann 1994) and that the flare is often generated at the 
bottom of the corona, this suggestion may not be reasonable. Alternatively, the pattern is 
such that the cool gas may be concentrated in a small region around the star (see Fig. 9(b)). 
In this case, the column density of this region should also be much higher than 2.5x10^° 
cm~^, since the best-flt Nh is the average value over the whole star. However, the flare 
may be located outside of that region with no absorption, so that the best-flt Nh value is 
consequently much lower. Because no rotational modulation has been detected, this cool gas 
region might be a belt running parallel to the star's equatorial plane. This pattern may also 
explain the Chandra LETGS observation (Ness et al. 2002). The spectra show the emission 
hne at 133 A, which is not expected with an Nh of 2.5 xlO^'^ cm~^. But it may appear if 
the absorbing material is concentrated in a belt around the equator and the line emission 
region is distributed all over the stellar surface. The belt might be reahzed by a disk of cool 
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gas originating from the stellar wind of the B type primary. 

5. Summciry 

We have analyzed the XMM-Newton observation of the eclipsing binary Algol which 
contains an X-ray dark B type star and an X-ray bright K type star. We fit the overall 
spectrum of the K star with a two-temperature thermal plasma model. The X-ray absorbing 
column density Nh of 2.5x10^° cm~^ obtained for the overall spectrum of the K star is 
significantly higher than the ISM column density measured through Na I line absorption 
to the optical continuum of the B star, which implies the existence of cool gas surrounding 
the X-ray emitting K star. However, the Nh value obtained for the fiare-only emission is 
consistent with the Na I line derived absorption, and thus leads us to suggest that cool gas 
is distributed in a belt like region which runs parallel to the star's equatorial plane. We 
cannot rule out that the elemental abundances of the spectrum's cool and hot components 
are identical, but the abundance of Fe is clearly different at the 4(7 level for the quiescent 
phase and significantly more than 5a for the flare phase, indicating that the hot and cool 
gases are not well mixed. We also observed a significant Fe abundance increase in the hot 
component by a factor of ~ 1.6 during the flare, which supports the idea of chromospheric 
evaporation. 

This project is supported by the Natural Science Foundation of China. This work is 
based on observations with XMM-Newton, an ESA science mission with instruments and 
contributions directly funded by ESA member states and the USA (NASA). 
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Table 1: Spectral properties of Algol (2T VMEKAL model) 
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- 12 - 



Table 2: Spectral properties of Algol (2T VMEKAL model). The top six rows are for the 
cool component and the lower six rows are for the hot component. 



No. 


C 


N 


O 


Ne 


Mg 


Si 


S 


Ar 


Fe 



0.11 
0.05 
0.09 
0.02 
0.12 
0.01 



+0.11 
0.11 

+0.10 
0.05 
+0.10 
-0.07 
+0.14 
0.02 
0.18 
-0.12 
+0.52 
0.01 



2.18 
1.93 
1.49 
1.65 
1.42 
2.35 



+0.01 
1.23 
+0.23 
0.71 
+0.13 
-0.62 
0.28 
0.87 
0.28 
-().,"')8 
+0.41 
2.35 



0.26 
0.23 
0.21 
0.22 
0.23 
0.28 



+0.01 

-0.08 
+0.03 
0.05 
+0.08 
-0.04 
+0.04 
0.11 
+0.03 
-0.03 
+ 0.10 
0.14 



1.04 
0.72 
0.73 



1 


< 0.37 


< 1.82 


< 0.14 


2 


< 0.29 


< 1.48 


< 0.13 


3 


< 0.39 


< 1.41 


< 0.42 


4 


< 0.48 


< 3.08 


< 0.48 


5 


< 0.79 


< 1.94 


< 0.12 


6 


< 1.00 


< 5.02 


< 0.41 



-0.42 
0.09 
0.17 
+ 0.13 
-0.16 

'-'•'^^-0.20 
Q 40+0.14 
^^•^^-0.20 

1 n5+"-^-^ 

-'-■'^'^-0.19 

CTT^ 

0.20 
+0.28 
0.46 
0.48 
-0.37 
0.78 
0.55 
0.82 
0.51 
+0.29 
0.36 



0.27 
0.26 
0.23 
0.22 
0.21 
0.23 



0.20 
0.50 
0.47 
0.86 
1.19 
0.70 



0.27: 
0.24 
0.37 
0.25 
0.31 
0.25 



0.03 
0.10 
0.08 
0.05 
0.03 
0.06 
0.08 
0.05 
0.09 
0.06 
0.13 
0.10 
+ U.lb 
0.09 
+0.14 
0.12 
0.22 
-0.10 
0.30 
0.25 
+0.27 
0.25 
+0.12 
0.16 



0.30 
0.24 
0.29 
0.23 
0.10 
0.33 



+0.06 
-0.13 
+0.09 
0.10 
0.10 
-0.10 
0.08 
0.09 
+0.13 
-0.10 
+0.25 
0.21 

•^■^^-0.04 

r, 90+0.13 

^■^*-0.05 

16+°" 
'-'■^^"-0.04 
-1+0.20 
-0.10 
c;+0.23 
-0.13 
K+O.OS 
-0.09 



< 0.08 

< 0.08 

< 0.06 

< 0.08 

< 0.15 

< 0.09 



< 0.64 

< 0.39 

< 1.08 

< 1.08 

< 1.12 

< 0.62 



0.30: 
0.45j 

o.ie"' 



< 0.22 

< 0.26 

< 0.18 

< 0.34 

< 0.25 

< 0.16 



< 0.,36 

< 0.30 

< 0.55 

< 0.27 

< 0.85 

< 0.29 



0.15 
0.15 
0.15 
0.13 
0.13 
0.16 



0.24 
0.28 
0.27 
0.42 
0.44 
0.23 



0.03 
).03 
).()1 
).02 
).()3 
).01 
).()2 
).01 
).()3 
).01 
).()4 
).02 

ruB- 

-0.00 
0.04 
0.02 
0.05 
-0.01 
0.05 
0.03 
0.06 
0.03 
+0.03 
0.03 



Tabic 3: Comparison of coronal elemental abundances from XMM-Newton (this paper), 
ASCA (Antunes et al. 1994) and Chandra LETGS (Schmitt & Ness 2004). 



Element 


This paper* 


This paper* 


Antunes et al. 


Antunes et al. 


Antunes et al. 


Schmitt & Ness 




Cool plasma 


Hot plasma 


Low State 


Medium State 


High State 




C 


'J-'J' -0.03 




n.a. 


n.a. 


n.a. 


< 0.04 


N 


1 «/!+"■ 1« 
^•^4_o.50 


'-'■^-'--0.00 


< 0.1 


< 0.1 


< 0.1 


2.00 





"-"•^^-0.03 




U.OU_Q Q4 


Q O1+0.03 
""-"^-O.OS 


24+0 03 
^•^*-0.03 


0.25 


Ne 


>J. 'O_o.05 


0.651°:?? 


„ 7f.+0.10 


, 99+0.08 
^■^''-O.OS 


1 ns+008 

l-U»_o.08 


0.95 


Mg 


(, 94+0.03 


n 90+0.09 
U-^O-0.07 


^■4^-0.06 
'-'■^'^-0.05 


n 64+0 05 
^■0^-0.05 


'^■^'-0.04 


0.50 


Si 


25+° °'5 
"•^3-0.05 


'J-^' -0.07 


65+0 05 
'-'■O3-0.04 


-0.03 


0.45 


S 




'^•'^^-0.03 


'-'■^^-0.07 


21+0 05 


o-o9l°:Sl 


n.a. 


Ar 


0.26lgg 




0.1010:10 


14+012 
"•1*-0.12 


0.1010:11 


n.a. 


Fe 


14+001 
<J.14_0.01 


„ 0, +0.02 
U.OJ-.Q.Ol 


0.3010:°! 


o-,+0.02 
-0.02 


o.32l°:°l 


0.20 



* Average of best-fit values of the 6 spectra 



Table 4: Best fit results for flare-only spectra 



No. 


t 

(s) 


Nh 


T 

(lO^K) 


z 

(%) 


EM 


flux 

/10-llerg-j 




DOF 


F-test 
probability 


1 
2 


6000 
5000 


< 0.60(2o-) 

< 1.58(2<7) 


63.6it:;^ 
eo.ilti 


55+0-12 
64+0-16 


r- OO+0.21 

0.00_g 

4 .;f;+0.25 
^.OD_g 23 


9.24 
7.94 


1.11 

1.02 


373 
264 


9.50 X 10-6 
5.71 X 10-3 
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Fig. 1.— EPIC PN and MOSl 0.3 to 10 kcV and RGSl 0.3 to 2.1 keV light curves of 
Algol as a function of time and phase respectively, with time bins of 100 s. The PN and 
MOSl curves have each been moved up by 1 cts/s and the RGSl up by 2 cts/s for clarity. 
The horizontal hne segment marked as "back" defines the time interval during which the 
background spectrum was taken for the analysis of the fiare-only spectra. 
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Fig. 2. — Time distribution of the spectral parameters temperature (top), emission measure 
(middle) and column density for the cool component (index 1, "A") and hot component 
(index 2, "*"). A'^^, which is unique for the two components, is marked with "O". 
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Fig. 3. — Time distribution of the abundances of C, N, O, Ne, Mg for the cool component 
("A") and hot component ("*"). 
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Fig. 4. — Time distribution of the abundances of Si, S, Ar and Fe for the cool component 
("A") and hot component ("*"). 
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Fig. 5. — Two-temperature thermal plasma model fits to spectrum ^4. Panel a is 
where the top curve is for the PN and the bottom curve is for MOSl and M0S2, 
for RGSl and RGS2. The EPIC spectra show the presence of an Fe 6.7 keV line. 



for EPIC, 
while b is 




Fe(cool) 



Fig. 6. — Fe abundance of the hot component versus the Fe abundance of the cool component 
obtained from spectrum #1 of the quiescent phase. The best-fit values are marked by a cross; 
the contours correspond to errors of la to 5a stepped by la. 
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Fe(cool) 

Fig. 7. — Fe abundance of the hot component versus the Fe abundance of the cool component 
obtained from spectrum #5 of the flare phase. The best-fit values are marked by a cross; 
the contours correspond to errors of la to 5a stepped by la. 
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Fig. 8. — The fiare-only spectrum #2. The upper curve is for the PN and the lower curve is 
for MOSl and M0S2. 
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Fig. 9. — Sketch of the cool gas distribution around the K star. As discussed in the text, we 
prefer pattern(b). 



